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We present a method to accelerate the dynamical evolution of multiqubit open system by employing dynami¬ 
cal decoupling pulses (DDPs) when the qubits are initially in W-type states. It is found that this speed-up evolu¬ 
tion can be achieved in both of the weak-coupling regime and the strong-coupling regime. The physical mech¬ 
anism behind the acceleration evolution is explained as the result of the joint action of the non-Markovianity 
of reservoirs and the excited-state population of qubits. It is shown that both of the non-Markovianity and the 
excited-state population can be controlled by DDPs to realize the quantum speed-up. 
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I. INTRODUCTION 

Quantum speed limit (QSL) is of particular interest and 
has attracted much attention in tremendous areas of quantum 
physics and quantum information, such as quantum communi¬ 
cation E! , quantum computation E-Hfl , quantum metrology 
16—11] and quantum optical control 1 1214201 . In these fields, 
quantum speed limit time (QSLT) is a key concept, defined as 
the minimum evolution time in which a system evolves from 
an initial state to a target state. QSLT denoted by tq$l de¬ 
termines the theoretical upper bound on the speed of dynam¬ 
ical evolution. And by first fixing the actual driving time r, 
the shorter QSLT tq$l, the greater the capacity for potential 
speedup will be. In this situation, tqsl = r means that the 
evolution is already along the fastest path and possesses no 
potential capacity for further acceleration [21]. 

The speed-up evolution of an open quantum system is 
of great significance to the robustness of quantum simula¬ 
tors and computers against decoherence 0221 12311 . Since all 
quantum systems are unavoidably coupled to their environ¬ 
ments, many works have been done on the derivation of the 
QSLT for open system dynamics E [24j-[30ll since that for 
isolated system dynamics ll3T1437ll . Generally, these bounds 
of QSLT could be divided into two categories: Mandelstam- 
Tamm type based on the Cauchy-Schwarz inequality and 
Margolus-Levitin type based on the von Neumann trace in¬ 
equality. Moreover, Deffner et al. demonstrated that non- 
Markovian effects can lead to a smaller QSLT and therefore 
speed up quantum evolution [27]. On the other hand, Xu et 
al. further proved that the mechanism for speedup is not only 
related to non-Markovianity but also to the population of ex¬ 
cited states of the quantum system under a given driving time 
fiftl. It is the competition between the non-Markovianity of 
reservoirs and the excited-state population of qubits that ul¬ 
timately determines the acceleration of quantum evolution in 
memory environments. A question may now arise how to de¬ 
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sign a feasible and effective mechanism to achieve a accel¬ 
eration process via modulating the non-Markovianity and the 
excited-state population . This question is of particular inter¬ 
est in the weak system-reservoir coupling regime, where the 
evolution possesses no potential capacity to accelerate with¬ 
out any operating on the system. However, there are few re¬ 
searches on this question [38, 39]. 

In this paper, we propose a new method to accelerate the dy¬ 
namical evolution of a multiqubit open system by employing 
dynamical decoupling pulses (DDPs). Generally, dynamical 
decoupling technology is used to overcome decoherence by 
averaging the unwanted interaction with environment to zero 
I40M5511 . Interestingly, we find that, DDPs can also accelerate 
the dynamical evolution, or in another word, reduce the time 
for performing an elementary logical operation. As a conse¬ 
quence, DDPs can greatly increase the number of operations 
within quantum coherence time of the system. This can help 
to improve the speed of quantum computers and communica¬ 
tion channels. Meanwhile, compared with a single qubit sys¬ 
tem, the multiqubit system has more abundant structures and 
would have more extensive applications in quantum informa¬ 
tion processing. This is the main motivation of our present 
work. 

In this paper, we consider a multiqubit open quantum sys¬ 
tem in which each qubit only interacts with its own reservoir, 
there is no direct interaction among qubits, and all reservoirs 
are independent with each other. For simplicity, we take the 
initial states of multiqubit system under our consideration to 
be W-type states. We find that the QSLT has nothing to do 
with the number of qubits and the initial-state parameters, but 
it can be modulated by the number of DDPs to realize the ac¬ 
celeration of quantum evolution. And this speed-up evolution 
assisted by DDPs can be achieved for both the weak system- 
reservoir coupling case and the strong system-reservoir cou¬ 
pling case. In order to reveal the essential reason of the ac¬ 
celeration, the roles of DDPs in the non-Markovianity and the 
excited-state population are studied in detail. 

The paper is organized as follows. In Sec. II, we present 
the physical model under our consideration and derive the dy¬ 
namics of multiqubit open system under DDPs. In Sec. Ill, 
we investigate the effect of DDPs on the evolution speed of 




2 


(a) 


JI JI JI JI JI 


0 T 2T3T 4T nT 





— |1> 
■i—10> 


Reservoir 


(b) 

Reservoir 1 Reservoir 2 





Reservoir N-l 



Reservoir N 


FIG. 1: (Color online) Schematic diagram of the physical system 
considered in this work, (a) A qubit interacts with a reservoir in 
the presence of periodic dynamical decoupling pulses (DDPs). The 
energy separation of the qubit is ato, and the periodic dynamical de¬ 
coupling is implemented by a train of instantaneous n pulses applied 
at the time T, 2T, 3 T, ■ ■ ■, and nT. (b) N pairs of independent qubit- 
reservoirs in which DDPs are applied to each qubit. 


the multiqubit system and indicate the speed-up mechanism 
of the dynamic evolution. Finally, we conclude our paper with 
discussions and remarks in the last section. 


II. PHYSICAL MODEL AND SOLUTION 


We consider N independent qubits interacting with their 
own reservoirs locally, and each qubit is controlled by a train 
of DDPs, as showed in Fig. |T] The dynamics of the A-qubit 
open system can be obtained by the use of time evolution of 
each qubit-reservoir pair with the following Hamiltonian |53i] 

H = (jJqct +ct- + y] (i)kfl\ak + gk (cr+cik + fljcr-) 
k k 

oo 

+ ^^6(t -nT)(r z , (1) 


where cjq is the energy difference between the excited state 11) 
and the ground state |0) of the qubit, (r L are the Pauli raising 
and lowering operators, and a' k {a.k) is the creation (annihila¬ 
tion) operator of the kth mode of the reservoir with the mode 
frequency w*. The coupling strength between the qubit and 
its reservoir is denoted by gk, which is further characterized 
by the effective Lorentzian spectral density (56] 


J(oj) 


l rod 2 
2n (u>o - a >) 2 + A 2 ’ 


( 2 ) 


where yo is the Markovian decay rate, A is the spectral width. 
Typically, yo < A /2 corresponds to the weak system-reservoir 
coupling regime, while yo > A/2 corresponds to the strong 
system-reservoir coupling regime. Without any operation on 
the system, the dynamics is Markovian in the weak-coupling 
regime, while the dynamics is non-Markovian in the strong¬ 
coupling regime. 


The last term in Eq. (Q} corresponds to the dynamical de¬ 
coupling control with a train of ideal instantaneous n pulses. 
T is the time interval between two consecutive pulses. Assum¬ 
ing that the width of each pulse is sufficiently short, it can be 
treated as 6 function. Under the action of each pulse, the qubit 
will rotate around the z axis by n. Then under the unitary 
transformation U (t) — T exp(-/ dt\ f & (6 ~ ) cr z ], 

Hamiltonian in Eq. ([[} gives rise to the following effective 
Hamiltonian 

Het T = OJ()(T + (T- + ^ (Oka\ak 
k 

+ ^(-1 ) n gk{o-+a k +a\o--}, (3) 

k 


where n = [t/T] is the number of applied pulses within the 
driving time t, denoted by the integer part of t/T. From the 
effective Hamiltonian in Eq. we can see that DDPs only 
changes the sign of the qubit-reservoir coupling strength pe¬ 
riodically. This is the physical principle why enough DDPs 
can eliminate the effective coupling between the qubit and its 
reservoir, and hence the decoherence of open system can be 
efficiently suppressed. 

Under the single-excitation approximation in a pair of 
qubit-reservoir, the reduced density matrix of the qubit at the 
time t can be exactly solved as 


Ps (0 = 


pii(0)*~ P\o(0)k, \ 

P* l0 ( 0)< 1 — Pii (0) k 2 j’ 


or in the form of Kraus operators 

Ps(t) = ^ K i Ps (0) Kj, 
1 = 1,2 


where the two Kraus operators are given by 


(4) 

(5) 


K x = *, |1> <1| +10X01, (6a) 

K 2 = VWiOXil. (6b) 

When t e [nT, (n + 1 )T), the general solution of k, reads [53] 


Kt 


e~ MI2 [2A n F h , + (1 + AA n ) F 2n ], A = 2y 0 , 
e~ A,/1 [A„ cosh (A n d) + B n sinh (A„ii)], A + 2yo, 


which is a real number. The coefficients on the right-hand side 
of Eq. 0 are given by 


_ A 2 T{gi’' + -p”) _ p'l+p'/ 4 F ln 

U 4 yj(AT) 2 + 4 ’ 2 " 2 A 2 T 

A n = a+m'/ + a_m" , B n — f3+m n + + /5_m n _ 


where we have introduced the following parameters 
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with d = yj A 2 - 2/lyo- £, = -^1 + sinh , and m ± = 
4 s inh(f)±£. 

From Eq. ©, we can find that the excited-state population 
of the qubit at the time t scaled by the excited-state population 
of the qubit at the initial time (i.e., the ratio of the excited-state 
population of the qubit between at the time t and at the initial 
time) can be written as the following simple expression 

Pt = *?, ( 10 ) 

where k, is given by Eq. ([T|i. It should be noted that what 
Eqs. (0 or Q describe is an amplitude damping channel 
model. From Eqs. |4} we can see that P, reflects not only 
the decay of excited-state population, but also the decay of 
quantum coherence of the qubit. 

In what follows, we will focus on the situation of A + 2y (l . 
In this case, P, becomes 

P, = e~ Al [A n cosh (A n d) + B n sinh (A n d)] 2 , (11) 


an operator A are defined as the eigenvalues of V AAA. For 
a Hermitian operator, its singular values are given by the ab¬ 
solute value of its eigenvalues. And the fidelity between an 
initial pure state po = l*AoX*Aol and a target state p T is defined 
as f T = (i^o|p T |i/f()). Moreover, the authors in Ref. further 
proved that the ML-type bound based on the operator norm 
(p - oo) provides the sharpest bound on the QSLT, so we will 
use the following bound to derive the QSLT of the A'-quhit 
open system, 


tqsl 


t| 1-/ t | 

.£ A < dt 


(14) 


Clearly, we can obtain the QSLT by calculating the fidelity f T 
and the operator norm of L t p t . 

For simplicity, here we assume that the A'-qubit system is 
initially prepared in the W-type states. 


|t/^o) — oq 1100... 0) -t- 1010... 0) + ... -t- |000... 1), 

(15) 


which indicates that P, can be modulated by the parameters of 
coupling spectrum yo and A, as well as the number of DDPs 
n. 

Because every pair of qubit-reservoir is independent and 
identical, the dynamics of the A'-qubit open system can be 
straightforwardly given by 

p, = X K=iV ? )]po [ <i^a)L (i2) 

where K P] (pj = 1,2) denotes the Kraus operator for the jth 
qubit in Eq. d6ji. Having obtained the density operator at any 
time, one can investigate the QSL of the A'-qubit open system 
as described in the following sections. 


III. QUANTUM SPEED-UP VIA DDPS 


In this section, we study the quantum speed-up of dynamic 
evolution of the A'-qubit open system by employing DDPs. 
More concretely, we study the effect of DDPs on the QSLT 
of this system. This is because that by first fixing the actual 
driving time, a decrease in the QSLT results in an increased 
ability to speed-up, whereas an increase in the QSLT means 
the decrease of the potential speed-up capacity. 

Now we begin with the definition of the QSLT in an open 
system. The QSLT is defined as the minimum time of a system 
evolution from an initial state po to a target state p T , which is 
governed by the time-dependent master equation L,p, = p t . 
With the help of the von Neumann trace inequality and the 
Cauchy-Schwarz inequality, the QSLT reads [27] 


tqsl 


11 ~ fr\ 

min ^E\,E T 2 ,El^ 


(13) 


where E T p = \ £ dt \\L,p,\\ p , and ||A|| p = (A 1 ’ + ■ ■ ■ + A?,) 1/p 
denotes the Schatten p norm with descending singular values 
of operator A, Ai > /L > • ■ ■ > A„. The singular values of 


where the normalization coefficients satisfy YiJ=i \ a j\ — 1- 
Submitting po = into Eq. ([T2l >. the reduce density 

matrix of the A'-qubit system at the time t reads 


Pt = Pi l*Ao>(<Aol + (1 - -PA |0.. .0)<0.. .0|, (16) 

where P, is given by Eq. (fill) , and the derivation of Eq. ( fl6t is 
shown in Appendix A. Hence, the fidelity between the initial 
W-type states po = |i^o)(i/'ol a nd the final state p T is found to 
be 


fr = Pt, (17) 

which means that the fidelity between the initial W-type states 
and the final state is equal to the excited-state population in 
the final state. 

Meanwhile, Eq. ([16]) makes L,p, = P,(|iAo> <^ol - 
|0... 0) (0 ... 0|), and (L,p, f = L,p,. Thus, the operator norm 
of L,p, can be obtained as the absolute value of its eigenvalues, 

di = |P,|. (18) 

Then, using Eq. (IT4l ). one can get the expression of the 
QSLT 


tqsl 


t(1 -Pt) 




(19) 


which reflects the speed of the dynamic evolution from the 
initial W-type states l^o) to the final state p T by an actual evo¬ 
lution time r. 

Interestingly, we note that the QSLT has nothing to do with 
the number of qubits N and the initial-state-parameter aj & 
This is largely because the characteristic of dynamics with 
the initial W-type state, as shown in Eq. dT6t . According to 
Eq. (ITTI) and Eq. (fl9l ). the QSLT could be modulated by the pa¬ 
rameters of coupling spectrum yo and A, as well as the number 
of DDPs n within the actual evolution time t. Thus, the num¬ 
ber of DDPs can be used as a control parameter to improve 
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FIG. 2: (Color online) The QSLT versus the number of DDPs n for a 
fixed actual driving time. The dash line and the solid line correspond 
to 7 = 0.2T (in the weak-coupling regime) and y = 5T (in the strong¬ 
coupling regime ), respectively. Other parameters are chosen as At = 
10 and A = 1. 



FIG. 3: (Color online) The excited-state population in the final state 
P T as a function of the number of DDPs n. The dash line and the 
solid line correspond to y = 0.2To (in the weak-coupling regime) 
and 7 = 5To (in the strong-coupling regime ), respectively. Other 
parameters are chosen as At = 10 and A = 1 . 


the speed of evolution. We shall numerically investigate the 
quantum speed-up in the following. 

In order to observe the speed-up effect induced by DDPs, 
let’s focus on the QSLT tqsl for a fixed actual evolution time 
of the system r. When the actual evolution time r is fixed, 
for the QSLT ratio tqsl/t < 1 the lower ratio tqsl/t , equiva¬ 
lently, the much shorter tqsl , the greater the capacity for po¬ 
tential speedup. And in the case that the actual evolution time 
equals the QSLT, i.e., tqsl/t = L the evolution is already 
along the fastest path and there is no likelihood of further ac¬ 
celeration. 

In Fig. [2] we plot the QSLT tqsl given by Eq. ( IT9t ver¬ 
sus the number of DDPs n for a fixed actual driving time r. 
From this figure, we can see that in the weak-coupling regime 
(yo = 0.2 A, the dash line), the QSLT monotonously decreases 
with the increase of the number of DDPs n. This clearly indi¬ 
cates the quantum speedup effect induced by DDPs. However, 
in the strong-coupling regime (70 = 5T, the solid line), the sit¬ 
uation is different. From Fig.[2]we can see that the QSLT in the 
strong-coupling regime first fluctuates in a small n region and 
then decreases continuously with the increase of n in the same 
as the weak-coupling regime. That is to say, there may ex¬ 
ist the speed-down effect in the strong-coupling regime when 
the number of DDPs is small. But anymore, only if the num- 




FIG. 4: (Color online) The non-Markovianity in a pair of qubit- 
reservoir for the fixed driving time r as a function of the number of 
DDPs n. The circle lines represent the non-Markovianity measure T. 
The dashed lines and solid lines correspond to the non-Markovianity 
components T e=0 and T 8=Jr/4 , respectively. Parameters are chosen as 
70 = 0 . 2/1 for the weak-coupling regime, y 0 = 5 A for the strong¬ 
coupling regime. At = 10. and A = 1. 


her of DDPs is large enough, DDPs could induce the shorter 
QSLT in both the weak-coupling regime and strong-coupling 
regime, so as to enhance the capacity for acceleration. This 
accelerating effect is especially obvious in the weak-coupling 
regime: in the absence of DDPs (i.e., n = 0), the actual evo¬ 
lution time achieves the QSLT, and there is no likelihood of 
acceleration. 

In order to understand the physical mechanism of the quan¬ 
tum speed-up assisted by DDPs, we turn to focus on the non- 
Markovianity of the reservoir within the driving time for a 
pair of qubit-reservoir [57] and the excited-state population 
of the qubits. Physically speaking, an increase of the distin- 
guishability of a pair of states during any time interval im¬ 
plies the emergence of the non-Markovianity. This can be 
interpreted as a flow of information from the reservoir back 
to the qubit system. In Appendix B, we have derived the 
non-Markovianity measure I characterized in terms of the 
amount of information exchanged between the qubit and its 
reservoir m. By use of the optimal initial-state pair of the 
qubit |«Ar ( 0 )) = cos(9|l) + sin 6 exp (iip) | 0 ) and |t /<2 ( 0 )) = 
sin 11) - cos flexp (/(/j) |0) with 9 e [0, n/2] , tp e [0, 27 t], the 
non-Markovianity measure in the time interval r can be ex¬ 
pressed as 

r = max{r e= o,r e= ^ /4 }, ( 20 ) 

where the two non-Markovianity components T fl= o and Te=n/A 
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are given by 

r e=0 = f P,dt , T fe7r/4 = I" — ^=dt. (21) 

Jo,p,>o Jo,P,>o 2 

Here the optimal initial-state pair has been numerically_proved 
to be either {| 0 >, | 1 » (0 = 0 ) or {|+>, |-» (0 = n/ 4) EMU. 
We can see that the non-Markovian dynamics appears only 
when P, > 0, and the QSLT in Eq. ( IT9l > can be reduced to 

TQSL = TtW' (22) 

1 ^ 1 -P T 

which indicates that the QSLT is determined by two quan¬ 
tities: the non-Markovianity Tg = o with the initial-state pair 
(| 0 ), 11 )} within the evolution time and the excited-state popu¬ 
lation in the final state P T in Eq. (flTt . Hence, the memory ef¬ 
fect of reservoir and the excited-state population become two 
essential factors for speeding up the dynamical evolution. 

As an illustration, we plot the influence of the number of 
DDPs n on the excited-state population in the final state P T 
of the qubit and the non-Markovianity measure of the reser¬ 
voir T in Fig. [3] and Fig. [4] respectively. On one hand, from 
Fig. [3] we can see that in the weak-coupling case (the dashed 
line), the value of P T of the qubit rises sharply under the ac¬ 
tion of DDPs in the small n regime, and it increases slowly 
in the large n regime. In particular, P T can approach to unity. 
These excited-state population changes well match with the 
QSLT changes as shown in Fig. [2] On the other hand, from 
Fig. Ufa) we can see that the non-Markovianity is very weak 
with the maximum T max - 0.15. This means the environment 
of the multiqubit open system under the action of DDPs is a 
weak non-Markovian reservoir. Hence, we can conclude that, 
although the non-Markovianity of reservoir is the necessary 
condition for speeding up l l2li 1281138 1. the excited-state popu¬ 
lation is the dominant mechanism of the quantum speed-up in 
the weak-coupling case. It is through controlling the excited- 
state population of the qubits that DDPs can accelerate the 
dynamic evolution of the quantum system under our consid¬ 
eration. 

In the strong-coupling case (i.e., yo = 5/1), the influence of 
the number of DDPs n on the excited-state population of the 
qubit P T and the non-Markovianity measure of the reservoir 
T is different from the weak-coupling case. The solid line in 
Fig.0reflects changes of the P T with respect to the number of 
DDPs while the circle line in Fig. 0b) describes changes of 
the non-Markovianity measure with respect to the number of 
DDPs. On one hand, from Fig. [3] we can see that the qubit al¬ 
most remains in the zero excited-state population for the small 
number of DDPs. This means that the excited-state popula¬ 
tion does not affect the speed of the dynamic evolution for 
the small number of DDPs. On the other hand, comparing 
Fig. [3 with Fig. 0b) for the strong-coupling case, we can find 
that the QSLT changing rule under the small n condition well 
matches with that of the non-Markovianity measure. So the 
non-Markovianity plays the decisive role on accelerating the 
dynamic evolution for the small n case in the strong-coupling 
regime. From Fig.[3]and Fig. 0 b). we can see that both of the 
P T and the non-Markovianity vary soundly with the number 



FIG. 5: (Color online) The dynamics of P, for different numbers 
of DDPs: ;i=0 (solid lines), n=5 (dashed lines), 77=10 (dot lines), 
71=20 (dot-dashed lines), (a) The case of the weak-coupling regime 
(yo = 0.2/1); (b) The case of the strong-coupling regime (y 0 = 5/1). 
2 = 1 . 


of DDPs in the large n regime of n > 10. Comparing Fig. [2] 
with Fig. [3] and Fig. 0b), for the strong-coupling case we 
can find that both of the excited-state population and the non- 
Markovianity play an important role for the quantum speedup, 
it is the competition between the excited-state population and 
the non-Markovianity that takes the responsibility to acceler¬ 
ate quantum evolution in the large n regime of DDPs. 

Finally, we discuss the control mechanism of the non- 
Markovianity of a reservoir by the use of DDPs. From Fig.0 
we can see that on the whole, the non-Markovianity increase 
firstly and then decrease later with increasing DDPs number n. 
This can be explained by the effect of DDPs on the dynamics 
of each qubit. As shown in Fig. 0 the excited-state popula¬ 
tion of the qubit P, fast oscillates as the number of DDPs n 
increases. Overall, DDPs not only makes the effective cou¬ 
pling between each qubit and its reservoir decrease, but also 
makes the information exchange between them more rapidly. 
It is the competition between these two effects that leads to 
the nonmonotonic behavior of T as a function of n. More¬ 
over, contrasting the two subfigures Fig. 0a) and Fig. 0b), 
we can see that the value of T in strong-coupling regime is 
greater than that in weak-coupling regime, while in the weak- 
coupling regime the dynamics will be Markovian (T = 0) in 
the absence of DDPs (n = 0). Besides, the optimal initial-state 
pairs are also different in two regimes for calculation of the 
measure of non-Markovianity. In the weak-coupling case, the 
optimal initial-state pair is proved to be {| 1 ), | 0 )} (T = Tg = o), 
while in the strong-coupling case, the optimal initial-state pair 
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changes from {|+>, |-» (r = r e=jr/4 ) to {|1>, |0» (r = r 9=0 ) as 
n is increased. But anyway, V s= o reflects the measure of non- 
Markovianity. 


2013CB921804 and the National Natural Science Foundation 
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IV. CONCLUSIONS 

In conclusion, we have studied the effects of DDPs on the 
QSLT of the A-qubit open system when the A-qubits are ini¬ 
tially in the W-type states. We have found that DDPs can 
be used to accelerate quantum evolution of multiqubit open 
system in both of the weak-coupling and the strong-coupling 
regimes. In the case of the weak-coupling between qubits and 
reservoirs, it has been shown that the capacity of the quan¬ 
tum acceleration monotonously increases with the number of 
DDPs. In the case of the strong-coupling between qubits and 
reservoirs, it is indicated that the quantum speed-up capac¬ 
ity is the same as that in the weak-coupling case when the 
number of DDPs is large enough. While when the number 
of DDPs is small, it has been found that quantum speed-down 
may happen in the strong-coupling regime between qubits and 
reservoirs. The essential physical mechanism for the speed¬ 
up evolution is that the non-Markovianity of reservoirs and 
the excited-state population of the qubits jointly determine 
the QSLT. Under the action of DDPs, the non-Markovianity 
of reservoirs and the excited-state population of the qubits 
vary so as to lead to the quantum speed-up of the multiqubit 
open system. The non-Markovianity of reservoir is the nec¬ 
essary condition for speeding up in both of the weak- and 
strong-coupling regimes. The excited-state population is the 
dominant mechanism of the quantum speed-up in the weak- 
coupling case while the non-Markovianity of reservoir is the 
dominant mechanism of the quantum speed-up in the strong¬ 
coupling case when the number of DDPs is small. Our find¬ 
ings in the present paper may lead to development of effective 
methods for the speedup evolution of multiqubit open sys¬ 
tems. These findings in the present paper could prove useful 
to accelerate quantum evolution of multiqubit open systems. 
In future work we will explore the possibility of realizing the 
quantum speed-up for a more broader class of initial states, 
interacted multiqubits, more complicate environments by the 
use of DDPs. 
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Appendix A: Derivation of Eq. JT77t 


In this appendix, we present a detailed derivation of quan¬ 
tum state of the /V-qubit system at an arbitrary time, i.e., the 
quantum state given by Eq. ( IT6l ). When the A'-qubit system is 
initially prepared in the W-type states 


It^o) — o'! 1100... 0) + q/ 2 |010... 0) + ... + o//v |000... 1) 



(Al) 


where |y) means that only yth qubit is in the state |1) and the 
other A-1 qubits are in the state |0). Then the density operator 
of the initial state can be written as 


Po = IfAoX^ol 

N 

= y h w oi + y aja * k ^ oi ( a2 > 

;=i j*k 


For the single qubit case, making use of Eq. © it is straight¬ 
forward to show 


2 

Ki 

o 

o 

Ki 

l 

1=1,2 




2 

Ki 

|0)<1| 


1=1,2 




2 


uxoi 


«=1,2 




2 

Ki 

nxn 

l 

i = l ,2 





10) <0|, (A3a) 

*|0><1|, (A3b) 

K/UX0I, (A3c) 

|1><1| + (1-*?)I0><0|. (A3d) 


For the A-qubit case, by the use of Eq. (IA3b we can 
obtain the actions of the Kraus operators on I /) (k\ ( j,k = 
1,2,3,-•• ,A) as follows 


Y K=,*x w] IX <*l [»£,<.«] = ly> (k\ , (j 4 k) (A4a) 

Y K=1^, I j> 01 Co] = Ij> 01 + (l - k}) |0...0)<0...0|. 


(A4b) 



Making use of Eq. (IA2L Eq. (1A4I) and Eq. ( 1 1 21 ). the re- expressed as 
duced density matrix of the A'-qubit system at time t can be 
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J 


Pt 


X ( 0 ] 2 \af I j) 01 + X I J) <*l 


H\4*2,-f*N 


j±k 


K=i<»] 


X \<*j\ 2 3 4 5 [*?I j) 01 + (l - IC?) 10 . • ■ 0 ) <0 ... 0 |] + k; X 1 j) <*l 


j*k 


Kjpo + (l-Kj) | 0 .., 0 )< 0 .,. 0 | 

p, l^o)<*Aol + (1 -Pt) 10 ... 0 )( 0 .,. 0 |, 

I- 


(A5) 


which exactly is the expression given by Eq. ©. 

Appendix B: The non-Markovianity measure 

The non-Markovianity measure we employ here is based 
on the amount of information exchanged between the open 
system and its reservoir, which is defined as [57] 

T = max I dtcr(t). (Bl) 

P 1,2(0) Jo-(r )>0 

where cr(t) denotes the rate of change of the trace distance 

cr(f) = -y-2) \pi(t),p 2 (t)\ , (B2) 

at 

and the trace distance is defined as 

£>\px{t),p 2 (i)\ = l -Tx\ Pl {t)-p 2 {t)\. (B3) 

Physically, <r(t) < 0 (T = 0) holds for all dynamical semi¬ 
groups and all time-dependent Markovian processes, while 
cr(f) > 0 (r > 0) for non-Markovian dynamics. It is should 
be noted that the maximum in Eq. (IB 11 1 is taken over all 
initial-state pairs of the system. Here we only consider the 
non-Markovianity in a pair of qubit-reservoir. Fortunately, 
it has been proved that, for the case of a qubit, the optimal 
initial-state pairs are always antipodal points on the Bloch 
sphere l58ll . So we can assume piXO) = |i/q ,2 (0)) (t/q ,2 (0)|, 

|i/q (0)) = cos 0 | 1 ) + sin# exp 00 ) |0) and \if/ 2 (0)) = sin 0 |l) - 


cos Pexp (iip) |0> with 9 e [0, tt/ 2] , e [0,27r], During the 
evolution time t, the optimal initial-state pair pi( 0 ) and p 2 ( 0 ) 
evolve to the final-state pair pi(f) and p 2 (t). According to 
Eq. (0), we can obtain the simple expression of the trace dis¬ 
tance between two final states at time t 

£>(Pi( 0 ,P 2 ( 0 ) = ^cos 2 (29)pi + siir(2 0)P„ (B4) 

What is more, as the Bloch sphere governed by Eq. © 
is rotation symmetrical with respect to the z axis, the most 
strongest oscillating direction is probably either the pole di¬ 
rection or any direction in the equatorial plane. In other words, 
the optimal initial-state pair should be either ]|0), 11)} (9 - 0) 
or {|+>, |-» (9 = ;r/4) with |±> = (|0> ± e'Xl>) / V2 jUS 
Substituting Eq. ( IB4b and Eq. ( IB2b into Eq. ( IB 11 1 we 
can obtain the simple expression of the measure of non- 
Markovianity within the driving time r 

r = max{r fl=0 ,r e=3r/4 }, (B5) 

where the two non-Markovianity parameters IXo and r H=njA 
are defined by 

r s= o = f P,dt, (B 6 a) 

Jo,p,>o 

X .,/4 = f -r—-j=df. (B 6 b) 

J o,A>o z yPt 
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